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ABSTRACT 

Aims. We study in this paper the ice composition in the envelope around intermediate-mass class I Young Stellar Objects (YSOs). 

Methods. We performed a spectroscopic survey toward five intermediate-mass class I YSOs located in the Southern Vela molecular cloud in 
the L (2.85^1.0 pm) and M (4.55-4.8 pm) bands at resolving powers A /AT = 600-800 up to 10,000, using the Infrared Spectrometer and Array 
Camera mounted on the Very Large Telescope-ANTU. Lower mass companion objects were observed simultaneously in both bands. 

Results. Solid PLO at 3 pm is detected in all sources, including the companion objects. CO ice at 4.67 pm is detected in a few main targets 
and one companion object. One object (LLN 19) shows little CO ice but strong gas-phase CO ro-vibrational lines in absorption. The CO ice 
profiles are different from source to source. The amount of water ice and CO ice trapped in a water-rich mantle may correlate with the flux ratio 
at 12 and 25 pm. The abundance of PLO-rich CO likely correlates with that of water ice. A weak feature at 3.54 pm attributed to solid CH 3 OH 
and a broad feature near 4.62 pm are observed toward LLN 17, but not toward the other sources. The derived abundances of solid CH 3 OH and 
OCN~ are -10 ± 2% and ~1 ± 0.2% of the HiO ice abundance respectively. The H 2 0 optical depths do not show an increase with envelope 
mass, nor do they show lower values for the companion objects compared with the main protostar. The line-of-sight CO ice abundance does 
not correlate with the source bolometric luminosity. 

Conclusions. Comparison of the solid CO profile toward LLN 17, which shows an extremely broad CO ice feature, and that of its lower mass 
companion at a few thousand AU, which exhibits a narrow profile, together with the detection of OCN~ toward LLN 17 provide direct evidences 
for local thermal processing of the ice. 

Key words. (Stars:) circumstellar matter - Astrochemistry - ISM: molecules 


1. Introduction 

Dust grains play an important role in the evolution of clouds 
from protostellar cores to circumstellar disks. Since dust grains 
are the main source of opacity, they control the thermal bal¬ 
ance of clouds. The surfaces of cold grains act as heat sink for 
highly exothermic reactions to occur (e.g., formation of Hi) or 
provide sites for atoms and molecules to freeze on. The freeze- 
out of molecules like CO is found to be important for regu¬ 
lating the gas phase chemistry of other species (e.g., ?; ?; ?). 
The frozen atoms and molecules accumulate on top of a re¬ 
fractory core (silicates, carbonaceous compounds) and form an 
icy mantle or react with other species to synthesize more com¬ 
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* Based on observations collected at the European Southern 
Observatory at La Silla and Paranal, Chile (ESO Programme 164.1- 
0605) 


plex molecules. The relative chemical composition of this ice 
mantle is well determined after three decades of studies us¬ 
ing both ground-based and space-borne telescopes. The core¬ 
mantle grain model is supported by spectropolarimetry stud¬ 
ies (e.g., ?). Solid H 2 O, CO, and CO 2 abound in most lines 
of sight where ices are detected (?). Sometimes, minor species 
such as CH 4 (-2%), HCOOH (~2%), OCN- (-0.2-1%) (e.g., 
?), and H 2 CO (3-6%) are found (e.g., ? ; ?; ?). By contrast, the 
presence of other minor constituents such as NH 3 is controver¬ 
sial; its abundance relative to water ice is likely less than - 10 % 
(e.g., ?; ?; ?; ?; ?). It should be emphasized that most of these 
studies refer to high mass young stellar objects. 

Solid methanol (CH 3 OH) is a particular case. It epitomizes 
the importance of molecular solids in the understanding of the 
gas phase chemistry. The radiative association of CH I and H 2 O 
is an inefficient gas-phase process that yields methanol abun¬ 
dances of - 1 (L 9 relative to H 2 while abundances of 10 7 - 
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10 6 have been found in hot cores (e.g., ?; ?). The prevalent 
view is that the high abundance of gas phase methanol come 
from the release of large amounts of frozen methanol formed 
on grain surfaces (?; ?; ?). Methanol ice abundance relative to 
water ice is found to vary from less than 3% w.r.t. water ice in 
quiescent regions up to 30% around massive protostars (?). If 
all the methanol ice is released in the gas phase, the abundance 
of methanol in the gas phase with respect to FL will amount 
to 3 x 10 - 7 - 3 x 10 6 assuming that the abundance of wa¬ 
ter ice with respect to FE is ~ 10 1 (e.g., ?). A similar situation 
is found for low-mass objects. ?) set a stringent limit of 5% 
of methanol with respect to water ice for sources located in 
the Taurus molecular cloud, while ?) found abundant methanol 
ice (14-25% of water ice) in 4 out of - 40 envelopes around 
protostars observed with the VLT. Possible formation routes of 
solid methanol are also disputed. The formation rate of solid 
methanol by hydrogenation of CO ice in the absence of energy 
input (i.e. hot atoms, UV or particle irradiation) measured in 
laboratory experiments remains controversial (e.g., ?; ?; ?). 

The advent of 8 m class telescopes equipped with large for¬ 
mat arrays opens up the opportunity to study large samples of 
low and intermediate-mass sources. We present here the first 
observations of molecular ice features in the L (2.8-4.1 pm) 
and M (4.5-5.1 pm) bands toward class I intermediate-mass 
young stellar objects (YSOs) located in the Vela molecular 
cloud complex. The spectra were obtained in the context of a 
large programme using the Infrared Spectrometer And Array 
Camera (ISAAC) mounted at the Very Large Telescope ANTU 
(VLT-ANTU) of the European Southern Observatory (ESO). 
Two major absorption features are observable with ground- 
based telescopes, along with some weak features. The first 
strong feature is centered around 3.01 pm (~ 3300 cnV 1 ) and 
is usually attributed to the stretching mode of solid FLO. The 
study of the solid-water profile has been used to better explore 
the ice structure in the water matrix (e.g., ?; ?; ?). The feature 
shows a broad excess absorption beyond 3.2 pm whose origin 
is still unclear, although scattering by the larger grains ( 0 . 1-1 
pm in radius) in the size distribution is the best candidate (?; 
?)• 

The other important feature is the solid-CO band at 4.67 pm 
(2140 cm" 1 ), whose profile is sensitive to the shape and size of 
the grains as well as the ice composition and temperature and 
is therefore a diagnostic of the evolutionary state and thermal 
history of ices (?). As soon as the grain is warmed to ^ 12-15 K 
by the luminosity of the object, CO molecules can diffuse into 
the ice and form new bonds, changing the morphology of the 
mixture and thus the profile at 4.67 pm, or they can sublime 
back to the gas phase (?; ?; ?). The mobility of CO and its 
high abundance in cold icy mantles also explain why it is a key 
species for surface reactions leading to polyatomic molecules 
such as C0 2 and CH 3 OH (?; ?; ?). 

The VLT-ISAAC observations are used to constrain the 
physical and chemical conditions in the envelopes of a sample 
of intermediate-mass stars. The paper is organized as follows. 
We present the objects of our sample in Sect. 2 and the obser¬ 
vations and data reduction procedures in Sect. 3. The results on 
the water ice band are described in Sect. 4.1. Evidence for solid 
CH 3 OH is presented in Sect. 4.2. The CO-ice band is presented 


in Sect. 4.3. A possible correlation between the abundance of 
CO embedded into a water ice matrix and the IRAS 12 pm/ 
25 pm color ratio is discussed in Sect. 5. Conclusion are pro¬ 
vided in Sect. 6 . These data complement the survey of CO and 
other species for a sample of ~40 low mass YSO’s by ??) ob¬ 
tained in the same programme. 



Fig. 1. CO (1-0) contour map of the Vela Molecular Cloud complex 
(?). The complex is composed of 4 clouds named A. B, C and D in 
the map. The positions of the observed sources are indicated and the 
filled circles are the location of the 50 other embedded young stellar 
objects. Indicated as well are the locations of the Vela Pulsar and the 
HII region 259,-4. The map is an adaptation of Fig.l of ?). 


2. Protostars in the Vela molecular clouds 

We chose to observe 5 intermediate-mass protostars with 
bolometric luminosities L t ~ 100-1000 L 0 , corresponding to 
masses of 2-10 M 0 (?). The characteristics of the objects are 
summarized in Table 1. These objects are the brightest sources 
in the L- and M-band taken from the sample of Class I YSO 
candidates in the Vela molecular cloud (?), hereafter LLN. The 
location of the sources is shown in Fig. 1, together with that of 
other candidates not studied in this work. Four of these sources 
have been observed by ?). The Vela molecular cloud is located 
in the Galactic plane, and is reachable only from observatories 
in the Southern hemisphere. It lies in an RA range (~9h) which 
was not covered by the Infrared Space Observatory (ISO) satel¬ 
lite during the entire mission to avoid damage from the infrared 
radiation from the Earth. 

The complex hosts three supernova remnants (Vela SNR, 
Puppis A and G266.2 -1.2), the latter discovered by recent X- 
ray satellites (?). The spectral energy distributions (SEDs) of 
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the objects, red IRAS colors and a near-infrared slope ,v = 
d k)g(AF A )/cl log(,i) > 0 as defined by ?) are characteristic of 
YSOs. Near-infrared imaging has revealed that the chosen tar¬ 
gets are the brightest members of their respective proto-stellar 
clusters (?). The companion stars observed serendipitously (see 
Sect. 3) are fainter class I objects. The entire complex was 
mapped in the CO / = 1 —» 0 transition by ?). The complex 
is composed of at least four Giant Molecular Clouds (A, B, C 
and D), with individual masses exceeding 10 5 M 0 , where active 
star formation occurs. These dark clouds are relatively nearby 
(700 + 200 pc), therefore giving the possibility to study the ices 
in a typical galactic environment. The location of the different 
objects and of their companions coincides with a maximum in 
the CO / = 1 —> 0 integrated emission and with the highest 
visual extinction. Extinction studies toward background stars 
have shown little foreground absorption from diffuse clouds 
with a maximum of Ay-3 A (?). Assuming that the threshold 
extinction for detection of the water ice feature in this cloud is 
~3.3, the amount of water ice that resides in the foreground is 
negligible. The exact threshold value is unknown for Vela, but 
it could be much higher than 3.3 (?). 

3. Observations and data reduction 

The observations were performed using ISAAC mounted on 
VLT-ANTU in Paranal, Chile, in the L (2.85-4.1 /an) and 
M (4.5-4.8 yum) band in lanuary and November 2001 under 
mediocre conditions. ISAAC is a cryogenic mid-infrared (1-5 
//m) imaging facility and grating spectrometer. The instrument 
uses a Santa Barbara 1024 X 1024 pixel InSb array. The en¬ 
trance slit of the grating spectrometer was set to 0.6", which 
results in a resolving power of 600 in the L band and 800 in 
the M-band (low resolution mode). Since the average seeing 
at Paranal in the infrared was ~0.6", little flux from the ob¬ 
ject was lost. The slit was oriented such that the main target 
and a second nearby object, when present, could be observed 
simultaneously. An example of an acquisition image at /.-hand 
is shown in Fig. 2. The companion object is located between 
5 and 10" from the main source. The large detector combined 
with the 0.6" slit allowed the entire L- and M- band to be ob¬ 
tained in one exposure each. The integration times of 30-40 
min per object and by band were chosen such that the S /N ra¬ 
tio reaches 60 in the wavelength regions of good atmospheric 
transmission, allowing an analysis of the profile. Any other ab¬ 
sorption features with an optical depth of 0.05 (3cr) can be de¬ 
tected as well. In complement, three sources LLN 13 (IRAS 
08375^4109), LLN 17 (IRAS 08448-4343) and LLN 19 (IRAS 
0847CM-321) were observed at a resolving power of 10,000 
(0.3" slit) in the M- band. A summary of the observations is 
presented in Table 2. 

The objects were observed with a chopping along the slit 
with a throw of typically 15" together with nodding. This tech¬ 
nique removes the majority of the sky and background noise. 
After or before each science target, the spectrum of a standard 
star was obtained with the same setting (BS 3185 of spectral 
type F6II and BS 3842 of spectral type G8II). The difference in 
airmass between the target and standard star is always kept less 
than 0.015. The stars used for atmospheric correction were cho¬ 


sen to be late-type photometric standards. The targets have ac¬ 
curate published photometric measurements (see Table 2) and 
those values were used to estimate the flux of the main and 
companion objects. The companions are 3.5 to 18 times fainter 
than the main sources, which are the brightest objects in their 
respective proto-clusters. The absolute flux calibration is uncer¬ 
tain by 10%. The standard star spectra were not modeled and 
therefore hydrogen absorption lines in the spectrum were not 
removed. In consequence, apparent emission features at 3.04 
yum (HI Pfe 5-10), 3.30 yum (HI PL) 5-9), and 3.74 yum (HI Pfy 
5-8) in the data may be spurious features from the rationing by 
the standard star. 



Fig. 2. Acquisition image in the stellar field of LLN 17 (IRAS 08448- 
4343. The orientation of the image (1024 x 1024 pixels) is North- 
South and the pixel scale is 1 "/pixel and corresponds to the length of 
the arrow. The negative points are the result of the chopping. The main 
target is called the source IRS 17-57 in ?) and the companion object 
is source IRS 17-40. 

The spectra were reduced using in-house software writ¬ 
ten in the Interactive Data Language ( IDL ). The data reduc¬ 
tion procedure is standard for IR observations (bad pixel maps, 
flat-fielding, etc.). The distortions of the raw spectra were cor¬ 
rected. In the L-band, the wavelength calibration was obtained 
from arc lamp measurements. In the M-band, because of the 
scarcity of halogen lamp lines in this wavelength range, the 
spectra were calibrated by comparing the strong atmospheric 
absorption lines with high signal-to-noise ratio spectra of the 
atmosphere above Paranal measured by the ESO staff. The data 
points in the strongest telluric absorption features show the 
lowest signal-to-noise ratio and were therefore dropped. This 
is particularly the case for the telluric gaseous methane absorp¬ 
tion band around 3.25 yum. Error bars are omitted in the flux 
plots for clarity for the L-band spectra, but the atmospheric 
transmission is provided for information in Fig. 3. 



4 


W.-F. Thi et al.: VLT-ISAAC spectroscopy of YSOs in Vela 


Table 1 . Characteristics of the sources. 



Source 

(1) 

Lnir 

(%) 

(2) 

Lmir 

(%) 

(3) 

Liras 

(%) 

(4) 

^submm 

(%) 

(5) 

Lbol 

(L 0 ) 

(6) 

^env 

(M G ) 

(7) 

A v 

(Mag) 

(8) 

N » 

(pc- 2 ) 

LLN 8 

IRAS 08211^1158 

5.3 

10.5 

53.3 

31.0 

141 ® 




LLN 13 

IRAS 08375-4109 

5.3 

6.8 

72.8 

17.7 

960 ® 

2.6 <b) 

> 40® 

800 <w 

LLN 17 

IRAS 08448-4343 

1.3 

4.0 

60.8 

33.9 

3100 ® 

6.4® 

~ 40® 

3400® 

LLN 19 

IRAS 08470-4321 

4.8 

20.3 

65.8 

9.1 

1600 ® 

3.5® 

~ 40® 

2400® 

LLN 41 

IRAS 09017-4716 



80.7 

19.3 

470 




LLN 20® 

IRAS 08476^1306 

0.5 

3.1 

66.8 

29.7 

1600 ,b> 

2.2 (b) 


3900® 

LLN 33® 

IRAS 08576^1314 

2.6 

12.6 

80.6 

4.2 

91«.) 




LLN 39® 

IRAS 09014-4736 

1.3 

5.7 

71.1 

22.0 

807'“' 




LLN 47® 

IRAS 09094-4522 

3.4 

13.3 

74.4 

8.9 

21 ® 





Notes: The dots refer to values which are not available in the literature. 

Column l^f: The fractional luminosities are computed by ?). The indices refer to the wavelength range 1-5 pm (column 1), 5-12 pm (column 
2), 12-100 pm (column 3) and 0.1-1 mm (column 4) for NIR, MIR, IRAS and submm, respectively. L submm are estimated by ?) when the 
millimeter observations do not exist. 

Column 5: Bolometric luminosity. 

Column 6: The envelope masses (M env ) are derived front 1 mm continuum flux by ?) for four of our sources. Envelope masses are only given 
for those sources for which millimeter observations have been carried out. 

Column 7: Extinction estimated by ?). 

Column 8: N * refers to the maximum stellar surface densities measured. 

Reference: "?), *?), c VLT-ISAAC data analyzed by ?). 


Table 2. Summary of the observations and infrared fluxes taken from the literature (The companion objects are located within 30" from the 
main source). 



Coordinate 

J 

H 

K 

L' 

M 

Standard 

Resolution 6 

Source 

RA (J1950) 

Dec (J1950) 

(Mag)" 

(Mag)" 

(Mag)" 

(Mag) 

(Mag)" 

star 


LLN 8 

08 21 07.4 

-41 58 13 

13.00 

10.12 

7.59 

4.79 

4.01 

BS 3185 

Low 

LLN 8b 

08 21 08.1 

-41 57 44 






BS 3185 

Low 

LLN 13 

08 37 30.8 

-41 09 29 

>14.80 

12.14 

8.95 

5.87 

5.09 

BS 3842 

Low, Med. 

LLN 13b 








BS 3842 

Low 

LLN 17 

08 44 48.8 

-43 43 26 

14.0 

11.60 

9.13 

5.77 

4.7 

BS 3185 

Low, Med. 

LLN 17b 








BS 3185 

Low 

LLN 19 

08 47 01.4 

-43 21 25 

14.7 

11.92 

8.88 

4.08 

2.73 

BS 3185 

Low, Med. 

LLN 41 

09 01 44.0 

-47 17 19 

>14.80 

11.51 

8.15 

4.5 

3.3 

BS 3185 

Low 

LLN 41b 

09 01 42.3 

-47 16 43 

14.1 

11.31 

9.41 

7.3 

6.9 

BS 3185 

Low 

LLN 20" 

08 47 39.4 

-43 06 08 

>14.90 

12.73 

10.80 

7.7 

6.2 


Med. 

LLN 33" 

08 57 36.8 

-43 14 35 

>16.20 

14.1 

11.16 

7.19 

5.7 


Med. 

LLN 39" 

09 01 26.4 

-47 36 34 

11.79 

10.16 

8.55 

5.64 

5.2 


Med. 

LLN 47" 

09 09 25.6 

-45 22 51 

11.82 

10.10 

8.90 

5.94 

4.8 


Med. 


Notes: 

" Photometry taken from ?). 

b Low resolution R =800 in the L and M band; Medium resolution R =10,000 in the M band only. 
c Observed by ?). 


4. Observational results 

4.1. Water ice 

The spectra of the main objects, converted into optical depth 
scale, are presented in Fig. 3, while the L-band spectra of the 
companion objects are displayed in Fig. 4 in flux scale. A local 


blackbody that fits the photometric data at 2.5 pm and 3.8 pm, 
measured by ?) was subtracted from the L-band spectra to ob¬ 
tain the optical depth scale. The low signal-to-noise ratio in the 
spectra of the companion objects prevents detailed analysis of 
the profiles, although column densities can be derived. 
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Wavelength (uml 

2.94 3.12 3.33 3.57 



Wavelength (um) 
3.12 3.33 



Wavelength (um] 

2.94 3.12 3.33 3.57 



Wavelength (um) 

2.94 3.12 3.33 3.57 



Wavelength (um] 

2.94 3.12 3.33 3.57 



Fig. 3. Optical depth profiles of the 3 /mi ice feature for the 5 objects (dotted line) compared to ice coated grains at different radii using 
laboratory spectra taken from (?). 


A broad absorption feature extending from 3600 to 2600 
cm -1 (2.85 to 3.8 /an) is observed in all spectra including the 
companion objects. This feature, centered around 3.1 gm, is at¬ 
tributed to the O-H stretching mode of FEO ice (e.g., ?). The 
detection of methanol-ice toward LLN 17 (IRAS 08448-4343) 


at 3.54 yum is discussed in Sect. 4.2. The band shape is compa¬ 
rable to that found in quiescent clouds, for example in the line 
of sight of Elias 16 (?), a field star located behind the Taurus 
cloud, which traces unprocessed ices. Table 3 shows the optical 
depths of the water band for the main and companion objects. 
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Wavelength (urn) Wavelength (urn) 

2.94 3.12 3.33 3.57 2.94 3.12 3.33 3.57 



3600 3400 3200 3000 2800 2600 3600 3400 3200 3000 2800 2600 

Wavenumber (cm 1 ) Wavenumber (cm 1 ) 


Wavelength (urn) Wavelength (urn) 

2.94 3.12 3.33 3.57 2.94 3.12 3.33 3.57 



IRAS 41b 


IRAS 8b 

i- 

Flux (Jy) 

o 

o 


Flux (Jy) 

o 

o 


0.01 


0.01 



3600 3400 3200 3000 2800 2600 3600 3400 3200 3000 2800 2600 

Wavenumber (cm 1 ) Wavenumber (cm’) 

Fig. 4. L-band spectra of the companion objects. The water ice absorption feature is seen in all spectra. The spectra are dominated by the 
spurious features due to a bad telluric transmission correction. The dashed lines indicate the adopted continuum used to estimate the water-ice 
optical depth. 


The optical depths toward the companion objects are grouped 
around 1.1 + 0.3, whereas a larger scatter is seen for the main 
objects. In the remainder of this section, we focus our analysis 
on the main objects. 

An extended wing to the main feature is present in all our 
sources. The identification of the carrier which gives rise to 
this long wavelength wing is a long-standing problem (?; ?). 
To compare the strength of the wing relative to the main fea¬ 
ture, the L-band spectra in optical depth scale were normal¬ 
ized to unity at 3.01 yum. It appears that the spectra can be di¬ 
vided into two groups. The first group is composed of objects 
where the long-wavelength wing centered at 3.25 /urn is rela¬ 
tively weak compared to the main feature. Three objects fall 
into this group (LLN 13, LLN 17 and LLN 41). Objects in the 
second group show a stronger wing relative to the main feature. 
Two objects in our sample (LLN 8 and LLN 19) pertain to this 
group. The similarity in the band profile among objects within 
the same group is remarkable (Fig. 7). In both groups, the en¬ 
tire profile was fitted empirically with a linear combination 


of two Gaussians centered respectively at 3.01 (FWHM- 0.14 
yum) and 3.25/rm ( FWHM-QAA yum). 

A few objects exhibit an additional broad and shallow ab¬ 
sorption feature at 3.47 fim. The optical depth of this feature is 
reported in Table 3. The 3.47 fim feature is also found in higher 
mass protostars (??). The broad 3.47 yum feature is attributed to 
an hydrate (e.g. ammonia hydrate, ?), in which other contribu¬ 
tors such as sp 3 carbon can be present (?). The optical depth of 
the main feature at 3.01 yum is compared to that at 3.25 yum in 
Fig. 8. In addition, data of a selected number of other low- and 
high-mass protostars are included. A strong correlation can be 
found among objects in each group. The optical depth at 3.25 
yum is found to be 0.27 and 0.45 times that at 3.01 yum for the 
two groups, respectively. The tight correlations suggest that the 
carrier(s) of the extended wing are strongly related to water-ice. 
However, the correlations may be fortuitous since the number 
of objects in the sample is limited. 

It has long been realized that light scattering by large ice 
grains leads to additional extinction on the long wavelength 
wing of the water band (e.g., ?). Several models of water ice 
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Table 3. Ice column densities. The errors are at 3 <x level and are dominated by systematics. 



Source 

r(H 2 Oice) 

Mce(H 2 0) 

(10 18 cm -2 ) 

Tice(H 2 0) 

(K) 

McefOCN - ) 8 
(10 16 cm -2 ) 

Mce(CH 3 OH) 6 
(10 17 cm -2 ) 

LLN 8 

IRAS 08211-4158 

0.55 ±0.10 

0.80 ±0.15 

10-40 

<1.5 

<1.0 

LLN 8b 

IRAS 08211—4158b 

1.2 ±0.3 

1.7 ±0.4 

10-40 



LLN 13 

IRAS 08375^1109 

1.4 ±0.3 

2.2 ±0.4 

10-40 

<0.5 

<1.0 

LLN 13b 

IRAS 08375^4109b 

1.1 ±0.3 

1.7 ±0.4 

10-40 



LLN 17 

IRAS 08448-4343 

2.46 ± 0.50 

3.61 ± 0.70 

10-40 

4.3 ± 0.5 

2.5 ± 1.0 

LLN 17b 

IRAS 08448-4343b 

1.4 ±0.3 

2.2 ±0.2 

10-40 

<1.5 


LLN 19 

IRAS 08470-4321 

0.72 ±0.10 

1.05 ±0.2 

10-40 

<0.5 

<1.0 

LLN 41 

IRAS 09017-4716 

1.26 ±0.20 

1.85 ±0.4 

10-40 

<1.5 

<1.0 

LLN 41b 

IRAS 09017-4716b 

0.8 ±0.3 

1.2 ±0.4 

10-40 

<1.5 



Notes: 

" The 3 cr upper limits are derived assuming Av=25 cm 1 and A=1.0 x 10 -16 cm -1 molec -1 . Upper limits derived from low resolution spectra 
are higher than that from medium resolution because of the possible contamination by CO gas phase lines. 

b The 3 a upper limits are derived assuming Av=14 cnT 1 and A=2.8 x 10 -ls cm -1 molec -1 . No upper limits are given for the companion 
objects because the signal-to-noise ratios are too low to provide scientifically meaningful upper limits. 


Table 4. Ice abundances relative to water ice. 


Source 

Mce(H 2 0) 

(10 18 cm -2 ) 

JVice(OCN-)/ 

AWfUO) 

(%) 

W ice (CH 3 OH)/ 

A',c,(H 2 0) 

(%) 

N ice (CO)/ 

iV ice (H 2 0) 

(%) 

LLN 8 

0.80 ±0.15 

<0.2 

<12.5 

<22 

LLN 8b 

1.7 ±0.4 




LLN 13 

2.2 ±0.4 

<0.2 

<4.5 

55 

LLN 13b 

1.7 ±0.4 




LLN 17 

3.61 ±0.70 

1.7 

6.9 

15 

LLN 17b 

2.2 ±0.2 

<1.1 


17 

LLN 19 

1.05 ±0.2 

<0.7 

<9.5 

4 

LLN 41 

1.85 ±0.4 

<0.1 

<5.4 

22 

LLN 41b 

1.2 ±0.4 

<1.9 


<15 


have been developed in which part of the long-wavelength 
wing is attributed to scattering due to large grains (greater than 
0.2 ftm). Our model is similar to that previously used by ?). 
The silicate core radius is fixed at a constant value of 0.1 gm. 
The grain core is coated with a water ice mantle whose thick¬ 
ness is allowed to vary to match the observed spectra. To sim¬ 
plify the problem, the ice is assumed to have a single tempera¬ 
ture. The absorption and scattering cross sections are computed 
using a Mie scattering theory program for coated-spheres (?). 
The optical constants provided by ?) for water ice and by ?) 
for the silicate core were used. The computed spectra for to¬ 
tal grain radii of 0.3, 0.4 and 0.5 gm are shown overlaid on 
the astronomical spectra in Fig. 3. The peak wavelength of the 
scattering cross section is red-shifted compared to that of the 
absorption cross-section and can therefore account for part of 
the long-wavelength wing. The shape of the long-wavelength 
wing strongly depends on the actual grain radius. The maxi¬ 
mum radius derived from this model is 0.4-0.5 gm and the ice 
temperature is below 40 K. 

The column densities were estimated by integrating numer¬ 
ically the laboratory spectra over the water band from 2.8 to 


3.8 /urn using the band strength for pure FLO ice of A — 2.0 x 
10 -16 cm molec -1 at 10 K measured by ?): 

/ r v dv 

(1) 

where r v is the optical depth at wavenumber v (cm -1 ) and 
A (cm molec -1 ) is the integrated absorption cross section per 
molecule (band strength). The temperature and column densi¬ 
ties of the ice giving the best fits are summarized in Table 3. 
Uncertainties in experimental band strengths for water ice are 
much lower than that in the determination of the continuum 
around the water ice band. The inferred column densities varies 
from 0.8 to 3.6 X 10 18 cm -2 . 

4.2. Methanol ice 

Only LLN 17 (IRAS 08448-4343) exhibits a strong methanol 
feature at 3.54gm superimposed on the solid water absorption. 
The analysis of the methanol feature requires the subtraction 
of the extended red absorption of the water ice and of the 3.47 
/urn. An artificial subtraction of the water ice wing is needed 
because laboratory mixtures containing methanol ice do not 
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take scattering effects into account. The red-wing continuum 
is composed of two gaussians whose characteristics are given 
in Sect. 4.1. The 3.47 gm feature is modeled empirically by the 
sum of two gaussians as well ( A a = 3.47 yum, FWHM- 0.12 yum 
and /1„=3.44 pm, FWHM=0A8 /an). The remaining feature is 
compared with mixed ices of H 2 O and CH 3 OH with different 
relative abundances. The peak position and the profile of the 
solid methanol bands are known to vary with the amount of 
water in the mixture. The presence of water shifts the peak to 
higher frequency and narrows the methanol band. The best fit 
is obtained with the mixture fLOTTLOH^ 10:1. The different 
components of the fit and the total are shown in the upper limit 
of Fig. 5. Possible improvement of the fit with a two phase ice 
is tested. A two phase ice mantle composed of pure methanol 
and a mixture H 20 :CH;OH= 1 0 :1 in the proportion 30% and 
70% respectively improves marginally the fit (Fig. 6 ). 

Adopting the integrated absorption coefficient valid for 
pure methanol (A- 2.8 x 10 18 cm molec -1 , ?), the derived 
column density of methanol is (2.5 + 1.0) X 10 17 cm -2 (3cr), 
corresponding to a relative abundance compared to water ice 
of ^ 6.9 + 2%. The search toward the other objects was un¬ 
successful with an upper limit of t(CH 30 H)< 0.02 (3<x), cor¬ 
responding to limits on the CH 3 OH ice abundances of - 1 x 
10 17 cm -2 or a limit on the abundance of <4.5-12.5%. The 
methanol ice abundances and upper limits are summarized in 
Table 4. The relatively high upper limits (4.5-12.5%) stem 
from the lower water ice abundance in the line of sight of most 
objects (AT FLO) < 2.3 x 10 18 cm -2 ) compared to that of LLN 
17 (JV(H 2 0) = 3.6 x 10 18 cm -2 ). 

The detection of solid methanol toward the intermediate- 
mass protostar LLN 17 (IRAS 08448-4343) with a relative 
abundance with respect to the water ice of =* 6.9% and the 
non-detection (< 5-10%) toward the other objects can be com¬ 
pared with the variable abundances toward high-mass proto¬ 
stars (from <1 % up to ~ 30% in RAFGL 7009S and W33 A; 
?; ?; ?; ?) as well as low mass protostars (<5 up to 25%, ?). 
Some of the large variations in the methanol ice abundance oc¬ 
cur between objects of similar mass located within the same 
cluster for low-mass star-forming regions such as Serpens. The 
upper limits are relatively high because of the low water ice col¬ 
umn density found toward most objects. Deeper limits for the 
LLN sources are needed to determine whether similar CH 3 OH 
ice abundance variations also hold for the Vela star-forming re¬ 
gion. 

4.3. CO ice 

In the M- band atmospheric window (4.5-4 .8 pm), one can ex¬ 
pect to detect the solid CO absorption feature and/or the ro- 
vibrational band of gaseous CO. A strong feature centered 
around 4.67 pm (2141 cm -1 ) corresponding to CO-ice is de¬ 
tected in three out of the five sources, namely LLN 13 (IRAS 
08375^4109), LLN 17 (IRAS 08448-4343) and LLN 41 (IRAS 
09017-4716) and one of the companion objects LLN 17b 
(IRAS 08448-4343b). Their spectra are plotted in Fig. 9 (open 
circles). The gas-phase CO ro-vibrational lines can introduce 
error in the interpretation of the solid feature. In particular, the 


Wavelength (urn) 

2.941 3.125 3.333 3.571 




Fig. 5. The upper panel displays the water ice feature toward LLN 
17. The adopted 2 gaussian continuum, the 3.47 pm feature and the 
methanol ice:water ice = 1:10 mixture are overplotted. The lower 
panel shows a zoom around the methanol feature. In this plot, the 2 
gaussian continuum has been subtracted. 


P( 1) line of gas-phase CO lies at the center of the CO-ice band. 
Medium resolution spectra obtained toward three objects (LLN 
13, LLN 17 and LLN19) are presented in Fig. 10. CO ice is de¬ 
tected in the medium , but not in the low resolution spectrum 
(Fig. 11) toward LLN 19. The cause of the discrepancy is the 
strong CO gas phase absorption lines, which render the detec¬ 
tion of the weak CO ice feature problematic. 

In contrast to the water-ice profile, the FWHM of the solid- 
CO feature differs from source to source (4.6-12.2 cm -1 ). The 
peak positions lie in a narrow range (2139-2141 cm -1 ). The 
FWHM seen toward LLN 13 (IRAS 08375^1109), LLN 41 
(IRAS 09017-4716) and LLN 17b (IRAS 08448-4343b) are 
narrow (4.6-5 .8 cm -1 ). A narrow width is often observed to¬ 
ward field stars and young embedded low-mass young stellar 
objects (?; ?). The source IRAS 08448^1343 (LLN 17) shows 
one of the broadest solid CO features ever observed (Av=12.2 
cm -1 ). We fitted the spectra using the decomposition outlined 
in ?) who showed that every line-of-sight can be fitted by a lin¬ 
ear combination of a narrow middle component (me), a broad 
red component (re) and a blue component (be). The compo- 
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Wavelength (urn) 

3.448 3.508 3.571 3.636 



Fig. 6 . The methanol feature toward LLN 17 is fitted by a single 
methanol mixture (CH 3 OH:H2O=l:10) and by a 2 component mix¬ 
ture (70% CH 3 OH:H 2 O=l:10 + 30% pure CH 3 OH). The improve¬ 
ment with the 2 component mixture is marginal. 

nents are Lorentzian for the red component and Gaussian for 
the blue and middle components. The central wavenumber and 
FWHM of each component and their interpretation are given 
in ?). The phenomenological fits to the spectra allow an as- 
trophysical classification of the sources while more classical 
fits with laboratory data help in understanding the grain man¬ 
tle composition as well as the grain shape and size source 
by source. The best fits are found by a Simplex optimization 
method (e.g., ?). The measurement errors are set constant with 
wavenumber so that an uniform weight is given to all data 
points. 

The optical depths and derived columns densities are sum¬ 
marized in Table 5, while the abundances relative to H 2 0 ice 
are given in Table 4. The solid CO column densities were ob¬ 
tained by integrating numerically the three components. The 
band strength for pure CO ice was measured to be A — 1.1 x 
ICC 17 cm molec 1 at 14 K (?; ?; ?). The uncertainties in the 
measured A value are not taken into account but are smaller 
than that introduced by the continuum subtraction. Although 
the band strength shows a 13% variation when CO is mixed 
with other molecules and a 17% variation with increasing tem¬ 
perature (?; ?), the same value for A is adopted to compute the 
column density for the three components. The errors of ~ 20- 
30% reflect the noise in the observed spectra and the telluric 
features removal. Also presented in Table 5 are the CO-ice up¬ 
per limits for LLN 8 (IRAS 08211-4158). Most of the CO is 
located in the gas for the latter object (see Fig. 10). The anal¬ 
ysis of the gas phase CO lines is presented in a separate paper 
(Thi et al. in preparation). The column densities of the mid¬ 
dle and red component were estimated by integrating the each 
band individually. Again, a unique band strength of 1.1 x 10 17 
cm molec -1 is assumed. 

The line-of-sight toward LLN 13 (IRAS 08375-4109) 
shows a high middle/red CO ratio typical of CO ice seen toward 
low-mass protostars and field stars. LLN 41 (IRAS 09017- 
4716) is also well matched by a high middle/red CO ratio 
mixture. In contrast, LLN 17 (IRAS 08448-4343) exhibits a 




Fig. 7. Normalized optical depths of the Vela sources fitted by a com¬ 
bination of 2 Gaussians, one centered at 3.01 yum, the other one at 
3.25 pm. The upper panel shows the normalized optical depth for 3 
Vela sources. The spectra are overplotted one on each other. The sim¬ 
ilarity between the 3 sources is so strong that it is difficult to separate 
the sources in the plot. The lower panel shows the spectrum of 2 other 
Vela sources. The spectra in both panels are fitted by a linear combi¬ 
nation of two Gaussians. 

stronger red component than a middle one. The CO-ice profile 
and optical depth is similar to that found toward the high-mass 
star GL961E (?) although the FWHM of ~ 12.2 cm -1 is among 
the largest ever observed. The optical depths and column den¬ 
sities derived from the fits with the phenomenological com¬ 
ponents are shown in Table 5, together with complementary 
data obtained from ?). The total CO ice column densities using 
fits with laboratory data and phenomenological components are 
similar within the errors, although the individual components 
show larger variations. 

4.4. The 4.62 /urn" XCN" feature 

Other weaker bands can also provide constraints on the ener¬ 
getic processes in the vicinity of young protostars. In partic¬ 
ular, the so-called XCN band near 4.62 fi m, contained in the 
CO survey, has been cited as energetic processes tracer since 









10 


W.-F. Thi et al.: VLT-ISAAC spectroscopy of YSOs in Vela 


Table 5. CO-ice band optical depth and column density derived from the phenomenological decomposition. 


Source 

\rtotal 

iV CO 

(10 17 cm- 2 ) 

T bc 

‘co 

—me 

‘CO 

T fC 

‘co 

M bc 

iy co 

(10 17 cm- 2 ) 

Ajmc 

iy CO 

(10 17 cm- 2 ) 

N vc 

iy co 

(10 17 cm' 2 ) 

yyme/ 

iy co / 

M rc 

iy CO 

LLN8 

< 1.8 

<0.10 







LLN 13 

12.15± 1.71 

0.37±0.12 

1.90±0.12 

0.36±0.12 

1.10+0.39 

6.53±0.15 

4.52±1.62 

1.44 

LLN 17 

531 ±2.22 

0.25±0.15 

0.18±0.15 

0.32±0.15 

0.72±0.48 

0.61±0.19 

4.04±2.10 

0.15 

LLN 17b 

3.83±1.74 

<0.51 

1.13±0.51 

<0.51 

<1.5 

3.83±0.58 

<6.39 

>0.6 

LLN 19* 

< 1.8 

<0.10 







LLN 19 £ 

0.43±0.05 

0.02±0.005 

0.03±0.005 

0.025±0.005 

0.07±0.02 

0 .10±0.02 

0.25±0.05 

0.41 

LLN 41 

4.09±2.16 

<0.15 

0.48±0.15 

0.18±0.15 

0.14±0.15 

1.64±0.18 

2.27±2.04 

0.72 

LLN 41b 

< 1.8 

<0.10 







LLN 20 rf 

5.19±0.96 

0.12±0.07 

0.38±0.08 

0.29±0.04 

0.35±0.20 

1.29±0.27 

3.55±0.49 

0.36 

LLN 33 rf 

7.85±0.81 

0.21 ±0.05 

0.94±0.09 

0.33±0.03 

0.61±0.14 

3.20±0.31 

4.04±0.37 

0.79 

LLN 39'' 

0.44±0.16 

0.004±0.003 

0 .02±0.01 

0.03±0.01 

0.011±0.008 

0.068±0.034 

0.37±0.12 

0.18 

LLN 47 rf 

< 0.36 

<0.06 

<0.02 

<0.01 

<0.17 

<0.07 

<0.12 



"The error bars are 30" level. 
*Low resolution spectrum. 
^Medium resolution spectrum. 
rf Data from ?). 




Wavenumber (cm') Wavenumber (cm') 




Wavenumber (cm') Wavenumber (cm') 

Fig. 9. Best fit to the observed spectra using the three components phenomenological decomposition. The observed data points are plotted as 
open circles. The thick line is the sum of the three components after convolution by the profile of the ISAAC spectrometer (R — 800).. 
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Fig. 10. Medium resolution spectra (R — 10,000) of three objects. The gas phase CO absorption lines are clearly visible and dominate the 
spectrum ofLLN 19. 
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Fig. 8. Optical depth at 3.25 /jm versus optical depth at 3.07/rm. The 
straight lines can be fitted to the data. The origin of both lines is (0,0). 
The data for sources which are not located in the Vela cloud are taken 
from the ISO-archive. 


Fig. 11. Continuum subtracted VLT ISAAC M-band spectrum of 
LLN 19 around the CO ice feature at 2140 cnr 1 (4.675 /jm). The 
spectrum is in an optical depth scale. The thick (red) line is the sum of 
the three fitting components. 
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its discovery toward the massive protostar W33A (?). The fea¬ 
ture is particularly strong toward LLN 17 (see Fig. 9), and is 
widely atttibuted to vibrational sUetching modes of -CN groups 
in molecular coatings on dust grains. The best candidate is 
the OCN ion (e.g., ?), as first proposed by ?). Quantitatively, 
adopting an integrated band sUength A = (4-10) x 10 17 cm 
molecule -1 for OCN - (?), the feature corresponds to a column 
density of 4.3-11 x 10 16 cm -2 . Recent laboratory work sug¬ 
gests however thatA maybe closer to 1.3 x 10 16 cm molecule -1 
(?). The concentration relative to FRO using the latter value 
is ~1%. The derived column density is similar to that found 
around other YSOs (e.g., ?; ?). Upper limits are difficult to es¬ 
timate for the other objects because the wavelength range of 
the OCN - feature is dominated by gas phase CO lines in emis¬ 
sion or absorption. The upper limits are given in Table 3 and 
the relative abundances in Table 4. Noteworthy, OCN - is de¬ 
tected in the object that shows the broadest CO feature, hence 
the most processed, but has one of the lowest CO/FUO ratios 
in our sample. Further discussion on OCN - in LLN 17 is post¬ 
poned to Sect. 5.4. 

5. Discussion 

5.1. Companion objects 

Water ice is detected toward four companion objects. The col¬ 
umn densities are lower than those toward the primary object, 
except for LLN8b. The water ice may be located in the outer 
part of the extended envelope surrounding both the main object 
and the companion. CO ice has not been detected but the sam¬ 
ple is too small and the upper limit too high to allow further 
discussion, except for LLN 17, which is further discussed in 
Sect. 5.4. 

5.2. Effect of pores on the water ice spectrum 

The ice layers on interstellar dust grains may be best simu¬ 
lated by background deposition in the laboratory that forms 
a porous ice (e.g., ?). A porous structure has several advan¬ 
tages over the compact structure as a candidate for interstellar 
ices. A porous ice can retain a significant amount of molecules 
such as CO. Moreover, by increasing the temperature, the ad¬ 
sorbed molecules at the surface can migrate into the bulk where 
they are trapped, i.e., they remain in the mantle even if the 
grain temperature exceeds the evaporation temperature of the 
species. The annealing increases the diffusion of molecules into 
the pores. This point is further discussed in Sect. 5.3. 

The effects on the optical constants of a water ice matrix 
due to isolated (i.e. not connected with each other) inclusions 
or pores can be simulated using an effective medium theory. 
One of the formulations is the Maxwell-Garnett approximation 
(e.g., ? for a detailed description). The generalized Maxwell- 
Garnett formula with ellipsoidal pores is as follows (?): 

ft ice(l — £ice)(l ~ L) + 6\ c e [L + (1 — T)6j cc ] 

&ff = - (2) 

L + (\ -L)e 1C c-/(l -e KC )L 

where ej cc is the complex dielectric function of the ice, / is the 
volume fraction of the inclusion (here the vacuum) and L is the 


depolarization factor. This factor is equal to 1/3 for spherical 
pores and lies between 1/3 and 1 for prolate spheroidal cav¬ 
ities (long axis along the normal to the layer and to the elec¬ 
tric field). By setting L — 1/3, the classical Maxwell-Garnett 
formula is recovered. The depolarization factor characterizes 
the effects of the shape of the pores but does not model the 
differences in the chemical bonding between a compact and a 
porous ice. The approximation was used to study the influence 
of various porosities and the depolarization factor on the sim¬ 
ulated spectra. The role of the depolarization factor is studied 
by keeping all other parameters constant. The resulting spec¬ 
tra are displayed in Fig. 12. The simulations were performed 
with a porosity of 0.5 (i.e. 50% of the volume is vacuum). It 
is difficult to estimate the porosity of actual interstellar water 
ice mantle, but 0.5 is most likely the highest possible value. 
The generalized Maxwell-Garnett approximation is in theory 
valid for porosities lower than ~0.3. Simulations with chang¬ 
ing the porosity from 0.2 to 0.5 does not result, however, in 
large differences in the shape of the profile. On the other hand, 
a small variation in the depolarization factor (e.g. from 0.3 to 
0.5) does. In Fig. 12, it can be seen that a higher depolariza¬ 
tion moves the peak of the absorption to shorter wavelengths 
and increases the absorption at higher frequencies compared to 
that at lower frequencies. This effect may jeopardize the esti¬ 
mate of the ice temperature in Sect. 4.1. The possible ambiguity 
between a low-temperature and high-temperature ice with dif¬ 
ferent depolarization factors is illustrated in Fig. 12. The simu¬ 
lated shapes are very similar in the blue. The simulations sug¬ 
gest that water ice between 10 and 40 K can fit the observed 
profiles. However, the 120 K optical constant of porous ice is 
not able to mimic amorphous cold ice. Porous ice at 120 K 
lacks absorption in the red but extra extinction can be provided 
by scattering by large grains. It should be noted that the water 
ice feature at 3/rm alone does not provide sufficient constraints 
on the porosity of the ice. Laboratory spectra of water ice at 
various porosity are warranted to test the validity of the model¬ 
ing. 


It is known that some warm water ice can be hidden in the 
broad profile (?). Crystalline HjO ice was observed in several 
lines of sight, mostly in the ejected envelope of evolved stars 
(??) and in the circumstellar environment of the massive YSO 
BN object (?). It is known that the ice mantle on grains around 
evolved stars is formed by condensation of water molecules 
synthesized in the gas phase (?), whereas the water ice mantle 
of interstellar grains is likely created by grain-surface reactions 
(?). In conclusion, a significant amount of high-temperature 
(Tice =40-60 K) ice can be hidden in highly porous ice because 
of this effect, making the estimate of the water ice temperature 
from the direct fit to the spectrum unreliable, although there 
is no direct evidence of highly porous water ice in the ISM. 
Further theoretical and laboratory investigations are needed to 
constrain the degree of porosity in ice mantles. 
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Fig. 12. The effects of the depolarization factor L on the water ice 
feature at 3300 cm 1 (~ 3.1 pm). The effect of porosity for water ice 
at 120 K is shown (/=0 means compact ice while /=0.5 means that 
50% of the volumn fraction is vacuum. The porous 120 K water ice 
absorption coefficient is close to the compact 40 K ice. 

U) 

TO 

5.3. Evidence for thermal processing of water and CO J= 
ices i 

The sources chosen for this study were selected among the w 
brightest class I young stellar objects found in the Vela molec- ^ 

ular clouds by ?). The objects posses similar characteristics but g 

exhibit a large variations in water and CO ice column densities. x| 
It is therefore interesting to find possible relations between the 
ice properties and the source characteristics. The source char¬ 
acteristics of interest are the extinction Ay, the infrared lumi¬ 
nosity L ir , and colors such as K - L, IRAS 12 p m/25 gm flux 
ratio R(12/25) and 25 yum/60 pm flux ratio R(25/60). The ra¬ 
tio R{ 12/25) is a good estimate of the warm dust temperature 
(100 < T < 250 K), while the ratio R(25/60) is sensitive to 
cooler dust (50 < T < 100 K). Finally, the K -L color allows to 
probe hot dust (T > 250). The source characteristics were cho¬ 
sen because they are available for all the sources in our sample. 

5.3.1. Water ice 

The water ice column densities are plotted against the infrared 
luminosity in the upper left panel of Fig. 13. They vary from 
0.8 to 3.6 x 10 18 cm -2 and do not correlate with increasing 
infrared luminosity measured between 7 and 135 pm. 

The extinction can be estimated from the observed E(H-K) 
with Ay= (15.3 ± 0.6) x E(H - K) (?). The intrinsic H - K 
color of early B up to F stars is close to zero. Assuming that 
the K band flux is dominated by the extincted central source, 

Edl - K) - II - K. The estimated extinction is an upper limit 
since the H - K color may be lower if emission from hot cir- 
cumstellar dust at ~ 1500 K dominates the near-infrared con¬ 
tinuum. Alternatively, the optical depth at 9.7 pm (Ay=(18.5 + 

1.5) x T 9 . 7 , ?)), measured in the IRAS-LRS spectra, provides a 
lower limit on the extinction because of potential intrinsic sil¬ 
icate emission. Therefore, the combined use of the two meth¬ 
ods brackets the actual extinction. The different estimates of 



Fig. 13. Water ice column density N{ FFO) versus the infrared lumi¬ 
nosity Lir. The error bars are 2 a level. 



R(12/25) 


Fig. 14. Water ice column density AffUO) normalized to the extinc¬ 
tion Ay versus the ratio R( 12/25). The error bars are 2a level. The 
errors on the ratios /?(12/25) are ~ 5-10 %. A trend is seen between 
A(H 2 0)/Av and tf(12/25). 

the extinction Ay are provided in Table 6. ?) estimated lower 
limits of 20 magnitudes for most sources in our sample. The 
discrepancy between the two methods is about a factor 3. In or¬ 
der to better estimate Ay, the full Spectral Energy Distribution 
has been modelled together with the IRAS-LRS data between 9 
and 25 pm of three objects (LLN 13, LLN 17 and LLN 19) us¬ 
ing the one-dimensional public radiative transfer code DUSTY 
(?) adopting a bare silicate grain model. The three objects were 
chosen because they have measured 1.3 mm fluxes. The best 
fits to the SEDs are shown in Fig. 15 and the estimated Ay are 
given in Table 6. DUSTY fails to fit the details of the silicate 
feature of (LLN 13, LLN 17), but provides reasonably good fit 
of the silicate feature for LLN 19. The absorption in 10-25 pm 
region of LLN 13 and LLN 17 can be caused by water and CO 2 
ice in addition to silicate absorption. The near-infrared and the 
SED fitting method give relatively close values for Ay. Thus, 
we will adopt the simpler near-infrared method in the rest of 
the discussion. The El- and K -band magnitudes are included in 
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Table 6. Estimation of the visual extinction by the 9.7 /mi optical 
depth, the near-infrared excess and the SED fitting method. 


Source 


Method 


” 19 . 7 ” 

”£(// - £)” 

”SED fitting” 

LLN 8 

> 8 

<39 


LLN 13 

>48 

<59 

60 

LLN 17 

> 13 

<46 

40 

LLN 19 

>22 

<56 

45 

LLN 41 

> 10 

<40 



Table 1. The derived extinctions using E(H - K) are between 
20 and 50, consistent with the values derived from the silicate 
feature. 

Adopting the relation 1 Vh=( 1.6 x 10 21 ) x Ay cuT 2 , the de¬ 
rived water ice abundance is (1.3-5.9) x 10 5 and is lower than 
the mean value found in quiescent molecular clouds in Taurus 
(7 X 10~ 5 , ?) and Serpens (9 + 1 X 10~ 5 , ?). The lower water 
ice abundances found toward the YSO’s in Vela may be as¬ 
cribed to thermal processing of the dust grains by the central 
object, lower average gas density, or by external radiation from 
surrounding YSO’s. ?) show that variations in the gas-phase 
abundances in the envelopes of YSOs can be explained by a 
variable size of the region over which the molecules are frozen 
out; a similar situation may apply here. 

To further ascertain the effect of thermal heating, the wa¬ 
ter ice abundance ATFLOj/A v is plotted against the IRAS 
12 fi m/25 fim flux ratio R( 12/25) in Fig. 14. The IRAS 12 gm 
band filter is relatively narrow and therefore, the IRAS fluxes 
are not strongly affected by the silicate absorption band. The 
water ice abundance clearly decreases with increasing warm 
dust temperature. This is consistent with the work of ?), who 
were the first to show that the color temperature of the dust 
correlates with ice heating/crystallization through the CCF ice 
bands, which trace CCE ice mixed with FLO ice. The heating 
of the water ice mantle at T >100 K should also results in 
the formation of a significant amount of heated water ice (nar¬ 
rower profile). However, the water ice observed in our sources 
is mostly in an amorphous state (see Sect. 4.1), although the 
narrow feature of warm water ice may be broadened in porous 
ices (see Sect. 5.2). Alternatively, the inspection of the IRAS- 
LRS spectra toward LLN 13 and LLN 17 shows that the 12//m 
photometry can be biased by water ice absorption. Indeed, the 
presence of large amounts of water ice with respect to silicate 
will decrease the ratio R( 12/25) and explain the trend seen in 
Fig. 14. 

5.3.2. CO ice 

In this subsection, we added four objects located in the Vela 
cloud studied by ?). A closer look at the source characteris¬ 
tics and their possible relationships with CO ice abundance are 
given in Fig. 16. In this figure, the abundance of CO ice for 
each object is represented by a circle whose radius is propor¬ 
tional to the CO ice column density. The upper left panel of 


Fig. 16 shows the line-of-sight extinction A v and infrared lu¬ 
minosity L ir of the nine sources. 

The total column density of CO ice does not correlate with 
the infrared nor the bolometric luminosity of the sources. In 
particular two sources LLN 17 and LLN 19 show similar values 
for Ay and Lir but CO ice has only been detected toward LLN 
17. 

It may be possible to constrain the dust temperature range 
to which the CO ice abundance is sensitive by finding relations 
between CO ice abundances and colors. The K - L, R( 12/25), 
and A*!25/60) colors encompass dust with temperature ranging 
from 50 to more than 250 K. This range of dust temperature 
is well above the CO ice sublimation temperature, thus we do 
not expect any correlation apart, from perhaps, a weak one with 
the coolest temperature /A25/60). From the two color-color di¬ 
agrams plotted in Fig. 16 (the upper-right and lower-left pan¬ 
els), it appears that the CO ice abundances vary with the ra¬ 
tio R( 12/25) but not with K - L, which is expected, nor with 
R(25/60), which is more surprising. Before discussing this find¬ 
ing, it is important to test whether the CO ice abundance is re¬ 
lated to the extinction in the line-of-sight. Water ice column 
densities are known to vary linearly with the visual extinction 
after a certain threshold value (e.g., ?). In other words, water 
ice abundances do not change with Ay. For the CO ice in our 
sample, the abundances do not correlate with Ay as testified by 
the lower-right panel of Fig. 16. 

Pure CO ice sublimes at ~20 K, while CO trapped in wa¬ 
ter ice, often associated with the red component at 2136 cm -1 , 
evaporates with sudden phase changes of the water matrix at 
higher temperature (??). The amount of CO ice trapped with 
water should therefore show a similarly decreasing trend as 
water ice with increasing value of R( 12/25). The column den¬ 
sity of the red component normalized to the estimated extinc¬ 
tion V rc (CO)/Av is plotted in Fig. 17. V rc (CO)/Ay decreases 
with increasing value of R( 12/25) . Both /V rc (CO)/A v and 
/VdLOj/A v show the same trend with R( 12/25) and therefore 
correlate with each other (right panel of Fig. 18). This corre¬ 
lation is consistent with the simultaneous sublimation of water 
and CO molecules trapped in the water matrix. It is also clear 
that N(CO)/A\r and /VfFLOj/A v do not correlate (left panel of 
Fig. 18). If R( 12/25) is another way to express VdLOj/A v , 
then Fig. 17 and Fig. 18 express the same relationship between 
the CO red ice component and water ice. 

It is difficult to solely attribute the red component to trapped 
CO molecules because of the absence in space of the 2152 
cnY 1 feature, which is present in laboratory data of CO-H 2 O 
mixtures. However, the 2152 cm -1 absorption feature seen in 
laboratory spectra of CO-water mixtures disappears when the 
dust is heated above ~80 K (?). Other effects may contribute 
to the red component. One possibility is that part of the red 
component is due to scattering by the larger grains in the grain 
size distribution similar to the effect seen for the red wing in 
the water ice band (see Sect. 4.1 and ?). In this scenario, the 
CO and water ice may located in two layers, which ensures 
that the 2136 cuT 1 feature does not appear even at low tem¬ 
perature. The CO ice trapped in a water matrix and the pure 
CO ice in large grains may be located in two separate popu¬ 
lations in the line-of-sight. The CO red component appears to 
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Fig. 15. Fit to the SED of LLN 13, LLN 17 and LLN 19. The upper panels show the fit to the entire SED while the upper panels focus on 
the near to mid-infrared range. Also shown are the contribution by the extincted central object, the dust scattered light and the dust thermal 


emission. 
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R(25/60) R(12/25) 

Fig. 16. The upper left panel shows the infrared luminosity of the sources computed using the IRAS 4 bands fluxes versus an estimate of the 
extinction toward the objects. The radius of the circles at the position of each object is proportional to the observed abundance of CO ice. The 
upper right and lower left panels are two color-color diagrams. The infrared ratios R(12/25) and RC25/60) are ratios between the fluxes F(Jy) at 
12 and 25 /mi and between fluxes at 25 and 60 pm respectively. The lower right panel is a extinction versus R(12/25) diagram. The dashed-line 
represents a possible diving line between presence and absence of CO ice. 


be a good temperature indicator. The possible relationship be¬ 
tween /V rc (CO)/.4 v and K( 12/25 ) should be tested in other star¬ 
forming regions. Likewise, future observations of ’COi ice 
will allow to test the possible relation between heated 'COi 
ice and the fraction of hot over cold grains for intermediate- 
mass YSOs. Spitzer observations of CO 2 ice toward the low- 
mass protostar HH 46 have shown that most of the ices are 
located in the cold (T < 50 K) part of a circumstellar envelope, 
although the inner envelope should be warm enough to explain 
the presence of a double-peaked absorption feature for the CO 2 
bending mode (?). The best fit to the CO 2 bending mode fea¬ 
ture is obtained with an ice mixture CFLOHiFLOiCC^O.3:1:1 
at 155 K in the laboratory or ~ 75 K in space. This indicates 
that the inner FLO, CO 2 and CH 3 OH ices are warm. 

5.4. Is LLN 17 a peculiar intermediate mass YSO? 

Among the intermediate-mass YSOs observed in the Vela 
Molecular Cloud, LLN 17 is the only object where large quan¬ 
tities of FLO-rich CO ice, of methanol ice and of OCN are 


found; although it should be emphasized that the relative abun¬ 
dance of methanol ice toward LLN 17 is not extremely high 
(- 6.9 % w.r.t. FLO ice). The ice abundances seen toward 
LLN 17 may be attributed to its peculiar environment. LLN 17 
is the most luminous YSO in our sample, it is located in the 
densest stellar cluster (see Table 1) and shows signs of molec¬ 
ular outflows (??). In this section, we discuss whether these 
particularities can explain the ice abundances. 

First we discuss the different ways to form solid-CFLOH 
and their implications. The chemistry of solid-CFLOH is still 
subject to considerable discussion. The exact moment of its for¬ 
mation and the amount of subsequent energetic processing are 
currently unknown. Originally it was suggested that methanol 
ice formed through hydrogenation of accreted CO (?), but re¬ 
cent laboratory experiments suggest that this mechanism has a 
very low yield (e.g. ?; ?), although the results are not conclu¬ 
sive (e.g., ?; ?). Hydrogen addition to CO at low temperature 
most likely occurs by H-tunneling since this reaction has an 
energetic barrier (?). The original scheme allows an efficient 
synthesis of methanol ice in cold pre-stellar cores where CO 
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Fig. 17. Red component column density V rc normalized to the extinc¬ 
tion A v versus the ratio R(12/25). The error bars are 2<x level. The 
errors on the ratios R( 12/25) are ~ 5-10 %. A clear trend is seen be¬ 
tween V r c/Av and R(12/25). 


molecule is known to be highly depleted onto grains although 
the amount of atomic H, which reacts with CO to form CH 3 OH, 
is also small (?). If methanol ice is formed at an earlier phase 
of star formation, probably at the pre-stellar core stage without 
significant energetic processing afterwards apart from the pro¬ 
cessing by UV photons generated by cosmic-rays interacting 
with Hi and H, similar abundances should be found for mem¬ 
bers of the same cluster, which is not the case for objects in the 
Vela molecular cloud and in the Serpens cloud (?). Likewise, 
methanol ice is not present in most lines of sight in the Taurus 
cloud (?). Therefore, methanol ice is probably synthesized at 
a later stage of stellar collapse where energetic events (UV 
and/or warm-up in shocks created from the interaction of stel¬ 
lar wind/outflow with the surrounding envelope) from the cen¬ 
tral object can provide the necessary energy to overcome the 
reaction barrier. LLN 17 is the most massive YSO in our sam¬ 
ple and, hence, the most UV luminous. It would be tempting 
to attribute the CH 3 OH ice to grain surface reactions triggered 
by UV photons from the central object. However, if we take 
into account all YSOs where methanol ice has been found, it 
is also clear that the methanol ice abundance does not correlate 
solely with the UV flux from the central star and direct stel¬ 
lar UV processing is also not dominant (e.g., ?; ?). In addition 
to the UV from the central star, shocks can produce copious 
amount of UV photons and heat the gas and dust. ?) have shown 
that large amounts of water ice can form in post-shock regions. 
According to this scenario, the water ice mantle is formed by 
condensation of water molecules formed in the high temper¬ 
ature region of shocks, which subsequently condense rapidly 
onto the bare silicate grains in the cool post-shock region. In 
theory, the water ice formed by vapor deposition at T > 100 K 
should be in the crystalline form. As discussed previously, crys¬ 
talline water ice is rarely seen toward YSOs. One possibility 
to circumvent the problem is that high water deposition rate 
can result in amorphous instead of crystalline ice (?). Another 
possibility is that crystalline ice is amorphized by cosmic-ray 


bombardment or Ly a radiation (?). As realized by ??), current 
gas phase chemical networks do not include all possible high 
temperature formation paths of methanol. As a consequence, 
the amount of methanol formed in shocks in standard mod¬ 
els is too low to account for the abundance in the solid phase. 
Another possibility to form efficiently CH 3 OH ice is that hot 
atomic hydrogens created by dissociation of H 2 in shocks im¬ 
pinge onto the grain surface to react with CO to form HCO, 
H 2 CO and finally CH 3 OH. Chemical models of shocks that in¬ 
clude all formation paths of molecules at high temperature are 
warranted. 

The methanol ice found in LLN 17 is most likely embed¬ 
ded into a water matrix (Sect. 4.2), indicating that water and 
methanol ice may have formed simultaneously or that methanol 
molecules have migrated. Interestingly, a collimated bipolar H 2 
jet structure composed of bright knots of line emission has been 
detected toward LLN 17 while neither LLN 19 nor LLN 13 
show sign of outflow (?). The morphology of the LLN 17 jet 
points to the presence of episodic phenomena of mass ejec¬ 
tion typically observed in protostellar jets (?). The molecular 
outflow has been mapped by ?) in the 12 CO 1-0 transition. 
Important unknowns in the outflow/shock scenario are the de¬ 
gree of mixing of the processed and unprocessed grains, the 
precession and the periodicity of the jet phenomenon. In sum¬ 
mary, detailed models are needed to determine whether shock 
chemistry is a viable mechanism to synthesize large amounts 
of gas and solid phase CH 3 OH. Likewise, correlation between 
outflow and high methanol ice content should be further inves¬ 
tigated. 



Fig. 18. Total CO column density IV(CO) and red component column 
density A rc normalized to the extinction Ay versus N( HiOj/Av in 
the left and right panel respectively. A clear trend is seen between 
A(CO) rc /A v and A(H 2 0)/A v . 

It is also interesting to relate the stellar and circumstellar 
characteristics of LLN 17 to the presence of OCN~ in its en¬ 
velope. The OCN ion is believed to be formed in a solid state 
acid-base reaction between NH 3 and isocyanic acid (HNCO), 
efficiently mediated by thermal processing (????). The main 
difficulty in this interpretation is that isocyanic acid ice has 
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never been detected in space, although it has a relatively high 
gas-phase abundance in hot cores and in shocked regions (?). 
The high abundance in shocked regions is consistent with the 
proposed outflow/shock scenario. 

Theoretically, HNCO ice abundance up to ~3% with re¬ 
spect to water ice can also be attained by grain-surface re¬ 
actions (?). Alternatively, an abundance of ~1% for OCN 
can be easily reached by UV-photolysis of the initial mix¬ 
ture H 2 O/CH 3 OH/NH 3 =100:15:15 (?). Although the required 
amount of UV-dose (fluence) exceeds the value estimated in 
dense clouds, the large amount of processed CO and CH 3 OH 
ice found toward LLN 17 suggests that the UV field around 
LLN 17 may be enhanced compared to that provided by cos¬ 
mic rays induced only. The abundance of NH 3 in the envelope 
around LLN 17 is difficult to estimate owing to the poor signal- 
to-noise ratio of the spectrum around the ammonia feature at 
2.97 fjm. However, sufficient amounts of NH 3 may have existed 
in the ice mantle since the 3.47 /vm feature is relatively strong 
toward LLN 17, assuming that the 3.47 yum feature is caused by 
ammonia hydrate. Noteworthy, the detection of OCN is often 
concomittent with that of CH 3 OH ice (see Table 5 of ?) and the 
references therein). 

Another particularity of LLN 17 is that the water ice abun¬ 
dance with respect to H 2 is the highest (see Fig. 18) in our 
sample. Methanol ice and OCN~ are often detected in objects 
with high water ice abundance (???). 

Finally, the presence of water-rich CO ice in LLN 17 and 
its absence in the companion object LLN 17b as seen in Figure 
9 reinforces the idea that the processing occurs in the vicinity 
of the central object (i.e. with a few thousand AU). The M- 
band spectrum of LLN 17b is too low to give meaningful upper 
limit on the amount of OCN . The difference between LLN 
17 and its companion suggests that external heating and/or UV 
processing by the other stars in the cluster play a minor role 
(see last column of Table 5). 

Interestingly, the OCN and methanol abundances relative 
to water in LLN 17 are close to those found in the envelope 
around HH 46 (?). Likewise, the water-rich CO abundance 
dominates over the water-poor CO and water ice is relatively 
abundant (5.7 x 10 5 ). 

In summary, the simultaneous presence of large amount of 
water-rich CO ice, CH 3 OH ice , and OCN ice toward LLN 17 
and the non-detection in other objects are consistent with the 
idea that these species are formed by a combination of UV and 
thermal processing in the inner regions of the circumstellar en¬ 
velope. The UV radiation can be generated by the interaction 
between the outflows and envelope, by the central source, or by 
a combination of both mechanisms. Another possibility is that 
CH 3 OH is first synthesized in the gas-phase, then condenses 
onto grains simultaneously with H 2 0 in post-shocked regions. 

6. Conclusions 

We have obtained L- and M- band spectra of a sample of 
intermediate-mass protostars in the Vela molecular cloud with 
the VLT-ISAAC. This is the first significant sample of interme¬ 
diate mass protostars for which ice data are published. 


A broad absorption feature at —3.01 /urn is detected in all 
sources (main and companion objects). The features show an 
extended wing beyond 3.25 /mi, which can be reproduced in 
part by scattering by grains at radius 0.4-0.5 gm. The water ice 
feature is dominated by absorption from cold amorphous ice 
although the spectroscopic signature of warm water ice can be 
masked if the ice is porous. 

Methanol ice is only detected around the protostar LLN 17 
(IRAS 08448-4343). The derived abundance is 10 ± 2 % rel¬ 
ative to water ice. The upper limit on the methanol abundance 
toward the other sources is between 5 and 10% with respect to 
water ice. 

Solid CO is detected in four main objects and one compan¬ 
ion object. The profiles show a large variety of shapes. A strong 
variation of the total CO ice column density is found. We de¬ 
compose the CO ice feature into three components. The column 
of CO becomes significant (i.e. larger than 10 17 cm~ 2 ) only at 
R( 12/25) greater than ~0.3. The color R( 12/25) may trace the 
abundance of water ice with respect to silicate. There is no clear 
trend between the column density of pure CO ice, traced by the 
middle component, and other characteristics of the YSOs (Lboi- 
Ay, IRAS bands flux ratios). On the other hand, we find a pos¬ 
sible correlation between the ratio of the flux at 12 /urn and 
25 yum, R( 12/25), which is a measure of the warm dust temper¬ 
ature (100< T < 250 K), and the amount of CO ice trapped in a 
water rich ice mantle, traced by the red component. Likewise, 
the amount of CO ice trapped in a water and that of water may 
correlate. This possible correlation is consistent with the idea 
that the water ice and the CO embedded in it sublime simulta¬ 
neously. 

These possible correlations should be tested in other star¬ 
forming regions. However, CO ice trapped in water-rich ice 
cannot solely account for the large amount of CO in the red 
component seen toward YSOs in Vela. Other factors such as 
scattering by the larger grains in the size distribution can prob¬ 
ably contribute to the red component. 

A strong absorption feature centered at 4.62 yum is detected 
toward LLN 17 (IRAS 08448-4343). The feature is likely 
caused by OCN A The derived abundance relative to water is 
1 + 0.2 %. This feature is not detected in any other object in 
our sample. Together with the detection of methanol and the 
broad CO feature, the detection of OCN suggests that the ice 
has been thermally and/or UV processed in LLN 17. The pro¬ 
cessing generated by UV from the central object is not essential 
and perhaps shocks induced processing is at play. Further the¬ 
oretical investigations on the possibility to form large amounts 
of methanol ice together with water ice in post-shocked regions 
are needed. 

Observations of a larger sample of high signal-to-noise ra¬ 
tio spectra obtained with 8-10 meter class telescopes (VLT, 
Gemini, Keck) and with Spitzer at longer wavelengths of pro¬ 
tostars of varying luminosities combined with sophisticated 
laboratory experiments will improve our understanding of the 
nature of the ices and their role in the synthesis of complex 
molecules in the interstellar medium. 
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